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We consider the adiabatic charge transport through zero-dimensional mesoscopic sample (quan- 
tum dot) caused by two periodically changing external perturbations. Both the magnitude and the 
sign of the transmitted charge are extremely sensitive to the configuration of the dot and to the 
magnetic field. We find the correlation function characterizing the random value of this pumped 
charge for arbitrary strength of the perturbation. In contrast to previous theoretical and experimen- 
tal claims, the charge is found not to have any symmetry with respect to the inversion of magnetic 
field. At strong pumping perturbation, the variance of the charge (Q^) is found to be proportional 
to the length of the contour in parametric space. 

PACS numbers; 72.10.Bg, 73.23.-b, 05.45. -fb 



I. INTRODUCTION 

Adiabatic charge pumping occurs in a system sub- 
jected to a very slow periodic perturbation. Upon the 
completion of the cycle, the Hamiltonian of the system 
returns to its initial form; however, the finite charge can 
be transmitted through a cross-section of the system. If 
the Hamiltonian depends on only one parameter, which 
is a strictly periodic function of time t, the value of the 
charge transfer, Q is zero. This may be not the case for 
Hamiltonians, which depend on two or more parameters. 
The pumped charge may be finite if these parameters 
follow a closed curve in the parameter space, which en- 
compasses a finite area. Such an evolution may be also 
characterized by a multivalued variable (angle); in this 
case transmitted charge is proportional to the winding 
number for this variable. The value of Q is not univer- 
sal. ThoulessB showed that for certain one dimensional 
systems with a gap in the excitation spectrum in the 
thermodynamic limit the charge Q is quantized. Such 
quantized charge pumping could be of practical impor- 
tance as a standard of electric currenilj. The accuracy of 
charge quantization depends on how adiabatic the pro- 
cess is. 

The practical attempts of creating a quantized elec- 
tron pump, are based on the phenomenon of Coulomb 
blockadeSa. In this kind of devices, several single elec- 
tron transistors (SET) are connected in series to increase 
the accuracy of charge quantization. At least two SET's 
are necessary to obtain a non-zero charge transfer. 

Semiconductor based quantum dotsU are often used 
now as Coulomb blockade devices. The advantage of 
these devices is the possibility of changing independently 
the gate voltage (and thus the average electron number 



in the dot) and the conductance of the quantum point 
contacts (QPC's) separating the dot from the leads. By 
doing so, one can traverse from almost classical Coulomb 
blockade to a completely open dot where the effects of 
the charge (wiantization are diminished. It was shown 
theoreticallyQ, that even weak backscattering in one- 
channel QPC's leads to almost quantized value of the 
pumped charge in low temperature limit, T 0. 

With the further opening of the QPCs the Coulomb 
blockade type charging effects become negligibleO. In 
this case, the main contribution to the pumped current 
is associated jAvith the quantum interference within the 
quantum domilj. Mechanism related to the inelastic 
processes was first considered in Ref. ^ In open sys- 
tems, this mechanism, however, is not effective, as the 
inelastic scattering time can be much larger than the 
dwell time of the electron in the dot. In this situation, 
the main mechanism of pumping is not related to the in- 
elastic processes but rather to change in-tJiei phase factors 
of the corresponding scattering matrixll3lll. Predictions 
of Refs. were apparently in accord with the recent 

data of Switkes et. alt2 

However, both papergl3'y are devoted to the weak 
pumping regime, in which the DC current is bilinear 
in the pumping amplitudes. On the other hand, in the 
experimental^ the pumping was not weak. Our purpose 
here is to construct a theory of the quantum pumping of 
a finite amplitude. The ensemble average of the pump- 
ing charge (Q) is equal to zero. We demonstrate that 
the variance of the pumping charge, (Q^) increases as 
the square of the area in the parameter space as long as 
(Q^) <C 1 . With further increase of the area, (Q^) in- 
creases much slower (as the length of the contour in the 
parametric space). These results are obtained in Sec. ||. 
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Another important problem to clarify iS|-tiie sensitivity 
of Q to applied magnetic field. Zhou et. a/.E3, claimed the 
symmetry of the pumped charge, similar to the Onsager 
relation for the conductanceoEj. The direct consequence 
of such symmetry is that the variance of of pumped 
charge is larger in the presence of time reversal symmetry, 
i? = than in its absence, ((3^(5 oo)) < ^Xi'^iB = 
0)) similarly to the conductance fluctuationsE3. Both 
those statemepis are supported by the experiment of 
Switkes et. alH. 

In our opinion, the statement about such symmetry 
is not correct. We discuss this problem in detail in 



Sec. ni. By manipulations with exact S -matrix we show 



that contrary to the conductance, the pumped charge 
is not symmetric with respect to inversion of magnetic 
field: Q{B) ^ Q{—B). Moreover, we show that at large 
values of the magnetic field the Q{B) and Q(—B) are 
not even correlated {Q{B)Q{—B)) 0. The variance 
of the pumped charge is found to be independent of 
magnetic field in agreement with the results of Ref. |ll| 
for large number of channels, and in disagreement with 
Refs. 1^,^. 

Our findings arc summarized in Sec. IV-i-Jn the same 
section, we will discuss existing experimentll3. 



II. MESOSCOPIC FLUCTUATIONS OF 
PUMPING IN ZERO MAGNETIC FIELD 



A. General formalism. 



The sample, therefore, is completely characterized by 
its unitary Nch x Nch scattering matrix, S{E), connect- 
ing ingoing and outgoing waves, with energy E. For in- 
stance, the two terminal conductance of the system G 
can be found from Landauer-Buttiker formula as 

G-^JdE[^-l^)Tr{nSiE)fJHE)}, (2.1) 



where f{E) = 1/(1 + e^/-^) is the Fermi distribution 
function (energy E is measured from the Fermi level). 
Matrices f/(r) are projectors on the states of left (right) 
lead: 



Q/3 



Sap X 



1, 



1< a< N, 



0, Ni<a< N,h ' 



(2.2) 



Let us now assume that S matrix is changing due to 
two external parameters Xi^2 slowly varying with time. 
The pumped charge Q is given byllil 

- ^ j dE (y-^ j dXidX2n {E, X) , (2.3a) 



n{E,X) = Im Tr<^ 



dS{E,X) dS^E.X) 



dX. 



dXi 



(2.3b) 



Our analysis will be based on the genergiLexpression 
for the pumped charge derived by Brouwenld based on 
the approach of Ref. n6l Consider the sample connected 
to two leads, as in Fig?]^. Each lead is characterized by 
its transverse modes a, where 1 < a < A^; for the left 
lead and Ni + 1 < a < Ni + Nr = Nch for the right lead. 




FIG. 1. Schematic picture of the sample connected to the 
leads "L" and "R". 



where X = {Xi, X2), and denotes the integration 
within the area encompassed by contour A. It follows 
from the unitarity of the S- matrix (charge conservation) 
that I1{E,X) can also be presented as 



niE,x) 



-Im Tr < T, 



dS{E,X) dS\E,X) 



dXo 



dXi 



(2.3c) 



Equations ( |2.3| ) are quite general. They assume only 
the absence of inelastic processes within the dot. We 
demonstrate the equivalence between the equations ( |2.3| ) 
and approach of Ref. ^ in Appendix ^ 

Closing this subsection, we discuss an important point 
— the physical meaning of the adiabatic approximation 
necessary for Eqs. (2^) to be valid. Usually, defining a 
perturbation of a quantum system as an adiabatic one 
means that the frequency lo of this perturbation is much 
smaller than the energy of the lowest excitations in this 
system. If the system were closed as in Fig. ||, the charge 
distribution after each period of the perturbation would 
return to the original distribution, and therefore, the 
pumped charge would be exactly quantized in units of 
the electron charge. 
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scattering phase 0: 



FIG. 2. Schematic picture of the quantum dot from Fig. |l| 
in the closed geometry. Genuine adiabatic approximation cor- 
responds to u • L ^ 0, whereas Eqs. (2.3) assume uj ■ L oo. 



However, for an open system, the spectrum is contin- 
uous, and the aforementioned adiabaticity criterium can 
be never satisfied. In this case we the only conditionwe 
can impose is smallness of the frequency as compared 
with the temperature and the mean level spacing Si of 
the dot separated from the leads. This condition is not 
sufficient for the qunatization of the charge. Neverthe- 
less, the value of the pumped charge can be still expressed 
in the form of the adiabatic curvature (|2.S 



To understand better the relation between closed and 
open systems, let us consider the gedanken experiment, 
where the two leadss are connected to each other as in 
Fig. ||. Level spacing of the whole system Ssys is propor- 
tional to 1/L. We will show now, that the result (2^) can 
be decomposed into two contributionsB. The first contri- 
bution is quantized and is just a charge which is pumped 
in truly adiabatic situation lu <C Ssys- This charge is, 
therefore, attributed to the dissipationless current. The 
second contribution is associated with the creation of the 
real electron-hole pairs in the system and is not quantized. 
It takes place in such closed system only provided that 
to ^ Ssys- This contribution involves dissipative conduc- 
tance and it is due to the fact that the electronic system 
can not fully adjust itself to the time evolution of the 
Hamiltonian: the state of the system at a given time is 
not an eigenstate of the Hamiltonian. Such retardation 
of the electrons from the external field leads to the dis- 
sipation. (Debye losses in closed systems with iv ^ Si 
gives a good example of such a dissipation). 

To illustrate such decomposition, we consider a simple 
case where each lead contains only one channel and there 
is no electron-electron interaction. In this particular case 
the S - matrix is 2 x 2 matrix which can be parameter- 
ized in terms of the dimensionless conductance g and the 



S = 



.id 



-ie 



(2.4) 



In the presence of the pumping perturbation, both pa- 
rameters g and 9 are slow f unc tions of time t. With the 
help of Eq. (2^), equation ( |2.3| ) gives 



Q — euj 



dt\\ 



, dQ 



(2.5) 



[This result was obtained by different means in Ref. 0] 
One can see th at t he first term in brackets of the in- 
tegrand in Eq. (2.5) gives always the quantized contri- 
bution because 0(l/27ra;) = (0) -I- 27rn, with n being 
an integer. At the same time, the term proportional to 
g{i) violates this quantization, and in the case of perfect 
conductance g = \ cancels the adiabatic contribution 
completely. 



B. Zero-dimensional model for the quantum dots 

In what follows we will be interested in the statistics 
of the pumped charge in the ensemble of quantum dots. 
In order to determine these statistics, we need to spec- 
ify the model. Firstly, we assume that the size of the 
quantum dots L is so small, that the Thouless energy 
Et ^ It erg far exceeds other energy scales of the prob- 
lem, such as the dephasing or escape rates (here T^^g is 
the characteristic time for the classical particle to cover 
all of the available phase space) . In this limit one can use 
the Random Matrix Theory to study the conductance of 
the system, see Ref. ^ All corrections to the RMT are 
small as Nchlgdot, where gdot — Et/Si and i5i is the 
mean level spacing. Secondly, we consider the adiabatic 
pumping in the dot with the large number Nch of open 
channels. In this approximation the effect of Coulomb 
interaction among the electrons in the dot turns out to 
be smaU as l/A^^h (^^^ R-^^- 1) ^'^'^ '^^'^ neglected. The 
condition Nch ^ 1 alsa allows us to use conventional dia- 
grammatic techniqucEj to perform the ensemble average, 
and to consider only lowest moments of the distribution 
of Q. 

The Hamiltonian of the system can be represented as: 



H = Ht 



Hr. + H, 



LD, 



(2.6) 



where Hd is the Hamiltonian of electrons in the dot, 
mimicked by a Af x M matrix 



Hr 



M 

E 

n,m— 1 



(2.7) 



We assume the "thermodynamic" limit M ~ g^ot —> oo. 
In the absence of pumping perturbations Hnm can be 
considered as a random (since gdot ^ 1) matrix that be- 
longs to the ensemble of real symmetric matrices (orthog- 
onal ensemble). Let the perturbations be represented by 
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two given not necessarily random (compare with Ref. EG) 



(1 2) 

M X M symmetric matrices Vn,m , 



so that 



Hnm{X) — Hum + XV nm — 'Hnm + -'^1 Kim + ^sV^m' 

(2.8) 

According to RMT, the correlation function of the ma- 
trix elements of the unperturbed part of the Hamiltonian, 
Hd, can be written as 

The coupling between the dot and the leads is 

Hld = J2 + H.c.) , (2.10) 

a,n,k 

where correspond to the states of the dot, tpa{k) de- 
notes different electron states in the leads (momentum k 
labels continuous spectrum in each channel a). 

The spectrum of electrons in the leads near Fermi sur- 
face can be linearized. Thus, without losing the general- 
ity we can write Hl as 



(2.11) 



where vp — l/inv is the Fermi velocity and v is the 
density of states at the Fermi surface. 



The coupling constants Wna in Eq. (2.1£1 are defined 



in the case of the reflectionless contacts asl 



M5i 



1, if n = a < Nch, 
0, otherwise. 



(2.12) 



For the system described above the scattering matrix 
S has the form: 

S^p{E,X) = 1 - 2niyWi^Gf,,^{E,X)W^p, (2.13) 

and the retarded (advanced) Green function G'^^{ Gj^^) 
is to be determined from the equation 



E - H{X) ± iTTiyWwA G'^'^i^E, X) = i. (2.14) 



Here matrices H and W are comprised by their elements 
( ^ and ( |2.12| ) respectively The factor in Eq. ( |2.12| ) 
is chosen so that the ensemble average scattering matrix 
{Sap) of a dot with fully open channels is zero. More 
complicat ed st ructure of W can be always reduced to 
the form ( 2.12 ) by suitable rotations. We are going to 
consi der t he averages of different moments of 11 given Jay 
Eq. (2.3b) with the help of diagrammatic techniquelia. 
For the technical reasons it is more convenient to trans- 



form Eq. (2.3b) with the help of Eq. (2.3c) 



n(£;,X) = Im Tr |f 

, NrTl - NlTr 



dS{E,X) dS^E,X) 



dXi 



N, 



ch 



(2.15) 



Notice that the matrix f is traceless. This fact signif- 
icantly simplifies further manipulations. We substitute 
Eq. (|2.13D into Eq. (|2.15D and obtain 



x"=x-^=x' 



sl^Tr {AG«(i?, X'^)VG^{E, X^)} 



(2.16) 



where = [e]*"' is the antisymmetric tensor of the sec- 
ond rank, and the factor s = 2 takes into account the 
spin degeneracy. We have introduced matrices A and F: 



A 



Notice, that 



1, 1 < n < Nch 
0, N^h <n<M 

1 <n< Ni 

Ni<n< N^h 
Nch<n< M 



(2.17) 




TrA = 0, Tr AF 



0. 



(2.18) 



To evaluate correlators of the functions Eq. (2.16), we 
adopted a diagrammatic technique for the ensemble av- 
eraging. In the thermodynamic limit M — > oo this tech- 
nique is somewhat-similar to the one developed for bulk 
disordered metalsliS. Factor 1/M plays now the same 
role as the small parameter l/epTimp with ep and Timp 
being the Fermi energy and the elastic mean free time 
correspondingly. The rules for reading those diagrams 
are shown in Fig. 



a) 




n 


n' 


m 




b) 





E + ( TtvW^W 



= xy,, 



d) 




© 



n m n n m n m n 



e) 
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FIG. 3. Elements of the diagram technique: (a) bare elec- 
tron retarded Green function; (b) correlator of the matrix 
elements of the Hamiltonian; (c) pumping perturbation; (d) 
renormalized electron Green function; (e) self-energy S. The 
second term in the self-energy, which includes an intersection 
of the dashed lines, is smaller compared to the first term as 
1/M. 

The ensemble averaged retarded Green function is 
given by 



1 



8 = 



E - H(X)±iTE 
1 

E - H{X) ± its - ■ 
M6i 



(2.19) 



where averaging is pe rfor med over realizations of random 
matrix Ti. from Eq. (2.8), and energy £ is of the same 



order as the width of the band of the random matrix 
eigenvalues. Self energy S includes, as usual, all of the 
one-particle irreducible graphs. In the leading in 1/M 
approximation it is given by the sum of the rainbow di- 
agrams, Fig. IsId, 



(2.20) 



Let us now expand the ensemble averaged Green function 
up to the second order in parametric perturbation. The 
final results are not nec essar ily quadratic in perturbation 
strength, see e.g., Eq. (2.25). All higher order terms can 
be neglected provided that 



\Vjli\«M('-l 



n,m — 1,2,. 



,M. 



(2.21) 



Inequality ( ^.2l| ) is nothing but the cond ition o f the cen- 
tral limit theorem. Solving Eqs. (2.19) and (2.2C), we 
find that at E £ 



±- 



Nrh ± ie 



4M 



(2.22) 



£2 



{XV) 



and first order terms, since the higher order terms vanish 
in the "thermodynamic" limit M —^ oo. 

Let us now turn to the analysis of th e statistical prop- 
erties of the function T from Eq. (2.16). All the relevant 
averages will be expressed in terms of the certain prod- 
ucts of the Green functions - diffuson V and Cooperon 

V (e, X) = Tr(G^' (ei + €,X + Y)G^ (ei, ^)) (2.24a) 
C (e, X) = Tr(G^ (ei +e,X + Y)G^ (d, ^)^), (2.24b) 



where we express energy in dimensionless units (2.23). 

The leading at M — > cxd and N^h ^ 1 approximation 
for the diffuson is a series of ladder diagrams. As usual, 
summation of this series can be performed by solving the 
equations presented graphically on Fig. ^ The solution 
of this diagrammatic equation is 



^ ,'27r 



1 



ie + Nch + iZ -X + XCqX 



(2.25) 



The diffuson has the univer sal form Eq. {^^) under the 
condition of Eq. (2.21). The structure and the strength 



of the perturbation potential V from Eq. {2.i) is encoded 
in two parameters: vector Z and tensor Cq. In terms of 
the original Hamiltonian they are given by 



Z 



Co 



2tt 
27r2 



-Try 



(2.26a) 



1,2, (2.26b) 



We used t he fa ct that the matrix V is symmetric. Pa- 
rameters (^.26 ) are also related to the typical value of 
the level velocities, which characterizes the evolution of 
energy levels of the closed system ei,{X) under the ac- 
tion of an external perturbation X ■ V, see Ref. |2^ (our 
definition is different by a numerical factor): 



dX 



Co 



— ( / ^^^^ 



(2.27a) 

(MM ■ 



where F is defined in Eq. (2.17). Here we introduced 
the dimensionless energy measured in units of mean level 
spacing 



27r£; 



(2.23) 



We can also expand these Green functions up to the first 
order in e/M and Nch/M and restrict ourselves by zero 



Now we are in a state to evaluate the statistics of the 
pumped charge. It should be noted that the specifics 
of the system enter only through the parameters and 
[Cojy. Moreover, all other responses of the system (e.g. 
parametric dependence of the conductance of the dai^ 
are also universal functions of the same parameters^. 
Therefore, Zi and [Co]ij can be in principle determined 
from independent measurements. 
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R 



A 



R 



+ 



A 



FIG. 4. Difluson diagrams. Cooperon diagrams differ from 
the ones for diffuson only in the direction of the arrow of the 
retarded Green function line. 



The ensemble average of the pumped charge vanishes 
and the sign of Q is random. To evaluate the typical 
value of this charge we consider the average {QaQb) for 
two different contours A and B on the parameter planes. 
To accomplis h thi s task, we need to average the product 
of functions ( 2.16 ) 




(a) 



i HB ; =: ; 




R(A) 

I 

R(A) 



R(A) 



R(A) 



(b) 



(2.28) 



where function T is defined by Eq. (2.1 



Diagrammatic expre ssion for I is shown on Fig. |[ Due 
to the relations ( 2. IS ), there is no need to renormalize 
the vertex A by the dashed lines. For the same reason, 
vertices A and F can not appear in the same cell and, 
there fore, the Cooperon Eq. ( ^.24bD does not contribute 
to I ( 2.28| ). We demonstrate in Sec. [II, that this implies 
the difference in how the change of sign of mag notic i fie Jd 
affects the pumped charge and the conductanccEJllj'ESl. 



FIG. 5. (a) Diagramm atic r epresentation for the correla- 
tion function X from Eq. (2.2J ); (b) Hikami box. 



The analytic expression for the diagram. Fig. ^ is 



Nch V 



(2.29) 



where diffuson pr opag ator T) is given by Eq. (2.25). The 
factor B in Eq. ( 2.2£ ) is given by the set of diagrams 
Fig. Ijb, which is analogous to the Hikami box for the 
disordered systems. It equals to 



S = Z?{1)+S(2)+S(3). 



(2.30) 



87r4 



51 



6 



SI 



Substitution of the expression for ^(2) into Eq. ( |2^ ) 
gives no contributio n to X due to the relationship be- 
tween X and n, Eq. (2.16). Indeed, the part of X which 
depends on B^^"^ does not contain . The contribution 
proportional to B^^^ vanishes in a similar way. Substitut- 
ing 6(1) into Eq. ( |2.16| ), and the result into Eq. ( p.29p , 
we find 



Uiei,X + Y)Uie2,Y) ) = 



1 



SNiNr 

N, 



ch 



92 



dX^dX, 



[ei- €2 - Z ■ X) 



XCqX 



(2.31) 



Here we used the explicit form of the diffuson (2.25), in- 
troduced dimensionless energies ( 2.26 ) and took the spin 
degeneracy into account. 

Now, we are in a state to evaluate the correlation func- 
tion of charges pumped in a course of motion along con- 
tours A and B on the parameter plane. The result be- 
comes compact if we choose new variables: 



X = 



y = 




(2.32a) 
(2.32b) 
(2.32c) 



Let us discuss why x,y, z are natural dimensionless 
variables. Recall that we are dealing with open sys- 
tems. Electronic escape time can be estimated as Te^c 
{Nch6i)~^- AH the energy levels have a finite width 
7 ~ I/tesc ~ NchSi > Si. It means that even at T = 
the pumping current is determined by the energy strip 
with a finite width 7 > (5i. The number of the levels 
in such a strip is a random function of the point in the 
parameter space. The correlation length R — \R\ of this 
random function can be estimated from the equation 



RZ -|- Ri 



Co 



R, 



7 



(2.33) 



The first term in the left hand side of Eq. ( 2.33| ) de- 
scribes a homogeneous shift of the spectrum, while the 
seco nd te rm represents random parametric oscillationsE3. 
Eq. ( 2.33|) can be rewritten in terms of the new variables 
Eq. ( |2.32[ ) as 



rz 



1 



(2.34) 



where r = {Co/NchY^'^R- As a result, r ^ min(l,Z^i). 
It means that in terms of the new variables x,y,z 
the pumping is weak (i.e. bilinear) provided that r < 
min ( 1 , z ~ 1 ) . Finally, 



{Q 



AWB 



d 



dxidx2 / dyidy2 



(2.35) 



dx 

a;_ = X — y 



2ttT 
NchSi 



where denotes the integration within the area encom- 
passed by contours a and b, see Eq. (2.32), and dimen- 
sionless correlation function /C is given by 



/C(u, -y, w) — 



dh- 



{2uh + vf + {1 + wf 



(2.36) 



f{h) = 



(/icotanh/i — 1) 
(sinh h)^ 



Low temperature regime corresponds to u <C 
max(l, V, w). In this limit 



IC(u, V, w) — 



1 



TT^ V^ + {1 + luY 



(2.37a) 



while at high temperatures u S> max(l, v, w) 



IC{u, V, w) 



1 



3ttu \1 + w 



(2.37b) 



Therefore heating suppression of the mesoscopic fluctua- 
tions of Q is similar to that of the conductance fluctua- 
tions. 



Equations ( 2.35 ) - ( 2.36 ) are the main results of this 
section. They describe the correlation between the charge 
pumped due to the motion in the parameter space along 
the different contours at arbitrary temperature. Now, we 
are going to apply Eq. ( 2.35| ) to analyze the variance of 
the charge {Q\)- 



C. Weak Pumping. 



If the characteristic magnitude of the potentials is so 
small that x'^_ <^ \, x^ ■ z <^ max (1, T/(7Vc?i(5i)), then 
the system is in bilinear response regime discussed in 
Refs. [lO|,[lll In this case one can put ?; = 0, w = in 



23§a 



Eq. (2.36D after the differentiation. As a result 



,2 , _ e'NiNrSl 



2nT 

NchSi 



z^lC-, 



2ttT 

NchSi 



(2.38) 



where Sa is the area enclosed by the contour a, in the 
parameter space x = {xi,X2)- Fu nction s ICi^2{x) can be 
expressed through f{h) from Eq. (2.36) in the following 
way: 
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ICi{u) = — / 



dh- 



f{h) 



u < 1; 



32 

T^, U > 1 



(2.39a) 



and 



dh 



f{h) (1 - Uu^h^ 



u > 1. 



(2.39b) 



In terms of the original pumping strength X Eq. (2.3S 
acquires the form 



(Qa) 



e^NiNrSl 



(det [Co] ) 



2tiT 

Nchh 



det 


Co 




1 


)^ 





ZjK,2 



(2.40) 
2nT 



Nchh 



Note that the specifics of the system enter onl y thro ugh 
the vector Z and the tensor Cq defined in Eq. (|2.26| ). In 
high temperature regime /C2 <C /Ci. It means that the 
simultaneous shift of all levels ( determined by z) is not 
relevant for pumping. On the other hand, at T — > 0, this 
simultaneous shift of all levels may be important. 



D. Strong Pumping. 



Let us now turn to the discussion of the opposite limit, 
3> 1, where pumping is strong. In this regime it is 



more convenient to transform Eq. ( 2.35| ) to the contour 
integrals. Using Stokes theorem, we find 



X =x y, 



dxi (p dyiK. 



2nT 

Nchh 



(2.41) 



We notice from Eq. ( 2.36 ) that the kernel /C decreases 
rapidly at a;^ > 1. Since the characteristic scale of the 
field itself is large a;^ 3> 1, we can perform the integra- 
tion over X- locally along the direction of the contour 
dx. It gives 



e^NiNr 



(h dxC 

J a 



2'kT z X X 



NchSi 



X 



(2.42) 



The kernel C can be expressed in terms of K. from 
Eq. (^!36|) as 



C {u, v) 



dwJC{u, vw, w^). 



(2.43) 



If the av erage d level vel ocity is small, y ^ 1, we find 
from Eqs. ( |2.36| ) and Eqs. ( ^.42D 



,0^^,2„ NiN ( 2^T \ 



(2.44) 



where is the length of the contour a. Function Ci in 
Eq. (|2.44D is given by 



25/2 






TT 




1 + 






4 

TT ' 


U < 1 


4 

3u ' 


U > 1 



dhf{h) 



l+4u2/i2 + (l + 4u2/j2) 



3/2 



1/2 



(2.45) 



In terms of the original pumping strength X, the dimen- 
sionless length of the contour ( Eq. (2.44)) acquires 
the form 



1 



N, 



ch J A 



(2.46) 



It is important to emphasize that in the case of the 
strong perturbation, the pumped charge is determined by 
the length of the contour rather than by its area, and it is 
not sensitive to the contour shape (provided that the con- 
tour is smooth on the scale of the order of unity). It has 
to be contrasted with the naive expectation (Q^) oc Sa, 
which follows from independent addition of areas. 

If z is not small, the value of the pumped charge de- 
pends not only on the length of the contour but also on it s 
shape . At lo w tem peratures, we obtain from Eqs. (2.42), 
( M3|) , and ( ^.37aD 



{Q\) 



dx- 



1 



[4 -I- 2;2 sin^ (zx)] 



1/2^ 



(2.47) 



where zx is an angle between the vectors z and x. At 
high temperature, (Q^) does not depend on on the shape 
of th e contour and is determined by Eqs. ( 2.44 ) and 



III. MAGNETIC FIELD EFFECTS ON 
ADIABATIC PUMPING. 



This section is devoted to the effect of the magn etic 

we 



IIIB 



field on the pumped charge. In subsection 
present a general discussion of the symmetries with re- 
spect to the time inversion. We demonstrate that, unlike 
the conductance, the pumped charge does not possess 
such a symme try. T his general conclusion is illustrated 
in subsection III B by the model calculation of the sec- 
ond moment of the charge pumped through the quantum 
dot. 
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A. Symmetry with respect to the reversal of the 
magnetic field 

Let us now consider the pumping through the meso- 
scopic sample subjecte d to a magnetic field B. The gen- 
eral formalism of Sec. II A remains valid. One can in- 



fer from Eqs. (2^) that the sign of the pumped charge 
changes together with the direction of the contour in the 
parameter space 



(3.1) 



where ^) denote opposite direction of motion in the 
para meter space along the same contour. Indeed, curva- 
ture (2.31:), is a single valued function of its parameters 
Xi, X2. Therefore, the only effect of reversal of the con- 
tour direction is to change the sign of the directed area 
dXidX2 — > ~dXidX2 without changing the integration 
domain. This immediately yields identity ( |3.l| ). Note 
that Eq. ( [3.1] ) relates the charges at the same m agne tic 
field. It is also important to emphasize that Eq. ( |3.l| ) is 
valid for arbitrary strength of the pumping potential. It 
is not restricted by the bilinear response regime. Eq. (23 ) 
of Ref. Q^iB) = -Q^{-B) together with Eq. ( [sl ) 
yield 

Q{B) - Q{-B), 

where the pumping is performed along the same contour 
in the parameter space. We intend,|te prove that, unlike 
for the two-terminal conductance!^ tI^h, such symmetry is 
not valid. 

The exact (not-averaged) S - matrix of the |S¥Stem 
changes with reversal of the magnetic field, B, asE^: 



S{E,B) = [S{E,-B)Y 



Sym metry relation for the t wo t eaninal conductance 
( prl) follows directly from Eq. (U)Ea 



G{~B) cx Tr{f,5(-B)f^5t(-B)| = 

= Tr{f,5^(B)f,5*(B)} ^ Tv {tJ {B)tiS^ (B)} = 

^Tr{TiS{B)frS\B)} 

We omitted factors independent of the magnetic field in 
the intermediate steps and used the unitarity of the S - 
matrix. Therefore, the relation 

G{B) = G(-B) 



is exact. 

Now let us turn to the pumped charge. 
Eq. (O) into Eq. (|23b|), we find 



n (-B) = Im Tr <^ n 



dS'^{B) dS*{B) 



dXo 



dXi 



Substituting 



(3.3) 



Im Tr { Ti 



dS^B) dSjB) 
dXi ' dX2 



From Eq. (3.3) one obtains 



n(B)-n(-S) = ^e'-'Tr<|f, 



dS^B)_dSiB) 



dX,, 



dX, 



(3.4) 



Commutator in right-hand-side of Eq. (3.4) vanishes only 
if ^-matrix is symmetric. This is not the case in the 
presence of a finite magnetic field. Therefore, there is 
no fundamental symmetry guarding the relation 11 (B) — 
n (— i?) as it was in the case for the two-terminal conduc- 
tance, and the Eq. (23) of Ref. |l^ does not hold. One can 
argue that it is the choice of the direction of the contour 
in the parameter space that violates the T - invariance 
and thus, B ^ —B symmetry. 

The absence of the symmetry with respect to the re- 
versal of the magnetic field, suggests that the correlation 
function {Q{Bi)Q{B2)) depends on the difference B1—B2 
only and vanishes at |-Bi — -B2I ^ 00 (as it does in generic 
parametric statistics). Model calculation of the following 
subsection confirms this expectation. 



B. Model calculation of the second moment 



In order to include the magnetic field into our descrip- 
tion we h av e to lift the condition that matrix Timn from 
(3.2) Eqs. (U|)-(U) is symmetric: 



where Ti-mn is the random realization of antisymmetric 
M X M matrix and a is the parameter proportional to 
the magnetic field. The resulting correlation function of 
two Hamiltonians at different values of the magnetic field 
Bi^2 can be conveniently presented in a form similar to 
Eq. (Ejl) as 



(7Y„„i(Si)7Y*,„,(S2)) - ( ^ 



(3.5) 



[Sn 



{AM 



D C 

Quantities ' characterize the effect of the magnetic 
field on the wave- functions of the closed dot and can be 
estimated as 



D 



Qdot 



$1 - ^2 



N?, = Qdot 



$1 -I- $2 
^•0 



(3.6) 
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where gaot ^ 1 is the dimensionless conductance of 
the closed dot, $1(2) is the magnetic flux through the 
dot which corresponds to the magnetic field i?i(2): a-nd 

$0 = hc/e is the flux quantum. 

For the Green functions in Eq. ( 2.24 ) taken at di ffer- 
ent magnetic fields Bi,B2, the diffuson. Eq. (2.25), is 
modified as 



V 



1 



^ie + Nch + Nj^ +iZ ■ X + XCqX 



(3.7) 



Diagrammatic repre senta tion, Fig. ^, and the expres- 
sion for H ikami box ( 2.30| ) remain intact. Instead of 
Eq. ( 2.35 ) we obtain 



{Qa{Bi)Qb{B2)) 
2 



2 , dxidx2 / dyidy2 

^^ch Ja Jb 
D 



x-=x-y, 



(3.8) 



Here the variables x, y, z. are deter mine d by Eq. ( 2.32| ) 
and the function K, is given by Eq. ( 2.3(: ). 

Equation ( |3.8| ) is the main result of this section. It de- 
scribes the sensitivity of the pumped charge to the mag- 
netic field. One immediately realizes that the correlation 
function depends only on the difference of the magnetic 
fields in accord with the discussion in the previous sub- 
section. Moreover the variance of the charge does not 
depend on the magnetic field (for. bilinear response this 
result was obtained by BrouweiO). In terms of the di- 
agrams, the absence of the symmetry with respect to 
the magnetic field reversal is revealed in the fact that 
the Cooperon does not contribute to the second moment 
even in the orthogonal case, [B = {]). 

We conclude this section by the discussion of the 
asymptotics of the correlation function {Q a{Bi)Q a[B2)) 
at finite magnetic field. We start with the weak pu mping 
at low temperatures T <C SiNch- Instead of Eq. (2.3S) 
one obtains 



{Qa{Bi)Qa{B2)) 



(3.9) 



2Nh 



where is related to the magnetic fields by Eq. ( ^.6[ ). 
If the dimensionless pumping potential is large (x'^ 3> 
[Nj^ /Nch + 1)) and the average level velocity is small, 
Z = 0, we obtain instead of Eq. ( 2.44 ): 



{Q\) 



fmNr 



Nr.. 



3/2 



(3.10) 



At high magnetic field, {Q a{B)Q a{— B)) rapidly de- 
creases as B~^, and one can use Eq. ( pj| ) for the bilinear 
response. 



Finally, we discuss the variance of the pumped charge 
in the magnetic field in the limit of high temperatures 



assuming the average level velocity is small, we find 



(T > max(2iii^^fJ:^,X2)). For the weak pumping. 



{Q\) 



ttT 



N, 



ch 



ch 



(3.11) 

In the strong pumping regime th e hig h temperature 
asymptotic behavior instead of Eq. ( ^.44 ) is given by: 



f2NiNr\ ( N,h6i 



N. 



ch 



Nr.. 



1/2 



(3.12) 



IV. DISCUSSION AND CONCLUSIONS 

Our main results include dependence on the pumping 
strength, temperature and magnetic field. 

Dependence of the pumping strength — At the small 
pumping potential, wEffissentially reproduced the results 
for bilinear responseEJllJ, that (Q^) oc S\ with Sa area 
being the area enclo sed by the contour in the parametric 
space, see Eq. ( 2.38 ). This bilinear response regime, how- 
ever, is valid only as long as the pumped charge is smaller 
than unity. The regime of strong pumping is analyzed 
for the f irst time in the present paper, see Eqs. ( 2.42| ) 
— (2.47). This regime is hallmarked by the dependence 
(Q^) oc Ia, with £a being the length of the contour, which 
is substantially slower than naive expectation (Q^) oc Sa- 
This slow dependence was already observed in Ref. [l^ . 
We think that our conclusion about independence of the 
pumped charge variancjc on the shape of the contour de- 
serves a careful checkEj. 

Temperature dependence — Our results for the high 
temperature regime T ^ Nch^i indicate that the vari- 
ance of the charge (Q^) is inversly proportional to the 
temperature (Q^) oc 1/T. Experimentt3 demonstrates 
(Q^) oc 1/T^ in the high-temperature regime. This dis- 
crepancy was attributed to the presence of the temper- 
ature dependent dephasijig, ignored in our treatment. 
In the simplistic modelsE^tJ, the dcphasing is described 
by adding an extra factor N^ into the mass of diffuson 



and Cooperon (3.7). If such a questionable procedure is 
adopted, the effect of dcphasing would be describ ed by 
replacement 



Nu, in formulas of Sec. IIIB 



We can see from Eqs. (3.11) — ( |3.12 ) that the same A^ 
produces different temperature depencleaces for the dif- 
ferent regimes. Use of experimentallycJ known depen- 
dence A^ oc T, would produce the results (Q^) oc 1/T^ 
\/T^/'^ for weak and strong pumping re- 



oc 



and (g2 

spectively. We believe that the available experimental 
information is not sufRcient yet for making detailed com- 
parison with our theory. 
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Effec t of t he magnetic field — We have demonstrated 
in Sec. Ill A that there is no fundamental reason for the 



pumped current to be symmetric with respect to the 
magnetic field reversal, in a sharp contrast with the de- 
pendence of conductance on the magnetic field. The 
corr espon ding correlation functions were calculated in 
Sec. Ill B| . It is demonstrated there that {Q{B)Q{—B)) cx 
at large B. These conclusions contradict to Ref. [ij 
where the symmetry with respect to magnetic field re- 
versal was reported. We can not explain this symmetry 
within the framework of our theory. 
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APPENDIX A: 

We demonstrate the equivalence between the equations 
( p.3[ ) and the approach of Ref. |l^. The demonstration 
will be based on equation of motion for the Green func- 
tions (2.14) related to t he S - matrices by E q. (2.13). 
First, we substitute Eq. (2.13) into Eq. ( ^.15 ), and ob- 
tain 



n(£;, X) ^ 2iTT^iy'^e'^ 



92 



(Al) 



xTr ^^WtW^G^ [e, X^) WW^G"^ [e, X''^) } . 

We will show now that matrix G^WW^G^ can be simpli- 
fied significantly using the equations for the Green func- 
tions (|1|). We pre-multiply Eq. ( |2l^ ) for G^ by G^, 
we transpose Eq. ( ^.14 ) for G^ and post- multiply it by 
G"*. Subtracting the results, we find 

2i7ri/G^ (e,X^) WW^G^ (e^X"^ 



(A2) 



= G^rB,X^ -G'^(E,X 



G^ E, X' 



G« [h(x^) 



Substituting Eq. (A2) into (|Al]), we obtain 
UiE,X) = 

e»^ -|-Tr LiyWrW^G'' (E, X) ^^G^ {E, X 



dX, 



dX. 



(A3) 

.1. 



If we recall that nvWrW^ is an operator of the current 
from the dot through the left contact, we obtain Eq. (10) 
of Ref. |l^, which proves that the physical mechanisms 
considered in Ref. fO and O are identical. 



1 D.J. Thouless, Phys. Rev. B 27, 6083 (1983). 

2 Q. Niu, Phys. Rev. Lett. 64, 1812 (1990). 

^ D.V. Averin and K. K. Likharev, in Mesoscopic Phenom- 
ena in Solids, Eds. B.L. Altshuler, P.A. Lee and R.A. Webb 
(Elsevier, Amsterdam, 1991). 

^ M.H. Devoret, D. Esteve and C. Urbina, in Les Houches. 
Session LXI Mesoscopic Quantum Physics Eds. E. Akker- 
mans, G. Montambaux, J.-L. Pichard and J. Zin- Justin 
(Elsevier, Amsterdam, 1995). 

^ L.P. Kouwenhoven et. al., in Proceedings of the Advanced 
Study Institute on Mesoscopic Electron Transport, Eds. 
L. Sohn, L.P. Kouwenhoven and G. Schoen (Kluwer, Se- 
ries E, 1997). 

^ L.P. Kouwenhoven et. al. Phys. Rev. Lett. 67, 1626 (1991). 
I.L. Aleiner and A.V. Andreev, Phys. Rev. Lett. 81, 1286 
(1998). 

* Effects of the Coulomb blockade on the conductance of the 
open dot were studied by P.W. Brouwer and I.L. Aleiner, 
Phys. Rev. Lett., 82, 390 (1999). 

^ B. Spivak, F. Zhou, and M.T. Beal Monod, Phys. Rev. B, 
51, 13226 (1995). 

F. Zhou, B. Spivak, and B.L. Altshuler, Phys. Rev. Lett. 

82, 608 (1999). 
" P.W. Brouwer, Phys. Rev. B 58, R10135 (1998). 
" M. Switkes, C.M. Marcus, K. Campman, and A.G. Gos- 

sard. Science, 283, 1907 (1999). 

B. L. Altshuler and D.E. Khmelnitskii, JETP Letters, 42, 
359 (1985). 

" L. Onsager, Phys. Rev. 38, 2265 (1931). 

M. Biittiker, Phys. Rev. Lett. 57, 1761 (1986). 

M. Biittiker, H. Thomas, and A. Pretre, Z. Phys. B 94, 

133 (1994). 

C. W.J. Beenakker, Rev. Mod. Phys, 69, 731, (1997). 

A. A. Abrikosov, L.P. Gorkov, I.E. Dzyaloshinskii, Methods 
of Quantum Field Theory in Statistical Physics, (Prentice- 
Hall, Englewood Cliffs, NJ, 1963). 

B. L. Altshuler and A.G. Aronov, in Electron-Electron In- 
teractions in Disordered Systems, A.L.Efros and M.PoUak 
eds., (North-Holland, Amsterdam, 1985) 

B.D. Simons and B.L. Ahshuler, Phys. Rev. Lett., 70, 
4063, (1993). 

21 B.L. Ahshuler and B.L Shklovskii, Sov. Phys. JETP 64, 
127, (1986) 

B.L.Altshuler and B.Z. Spivak, JETP Lett. 42, 447 (1985). 
^■^ L.D. Landau and E.M. Lifshitz, Quanutm mechanics, 3rd 

edition, (Pergamon Press, 1977). 
■^■^ Detailed comparison with the data was performed recently 

by M. Vavilov (private communication). 

H.U. Baranger, P.A. Mello, Phys. Rev. B, 51, 4703, (1995); 

P.W. Brouwer and C.W.J. Beenakker, ibid., 51, 7739, 

(1995); ibid., 55, 4695, (1997); I.L. Aleiner and A.I. Larkin. 

ibid., 54, 14423 (1996); E. McCann and I.V. Lerner, ibid., 

57, 7219 (1998). 

A.G. Huibers et. al, Phys. Rev. Lett., 81, 200 (1998). 



11 



